Oscillations in finite Fermi systems 
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Abstract 

A semiclassical linear response theory based on the 
Vlasov equation is reviewed. The approach discussed 
here differs from the classical one of Vlasov and Lan- 
dau for the fact that the finite size of the system is 
explicitly taken into account. The non-trivial prob- 
lem of deciding which boundary conditions are more 
appropriate for the fluctuations of the phase-space 
density has been circumvented by studying solutions 
corresponding to different boundary conditions (fixed 
and moving surface). The fixed-surface theory has 
been applied to systems having both spherical and 
spheroidal equilibrium shapes. The moving-surface 
theory is related to the liquid-drop model of the nu- 
cleus and from it one can obtain a kinetic- theory de- 
scription of surface and compression modes in nuclei. 
Quantum corrections to the semiclassical theory are 
also briefly discussed. 



1 INTRODUCTION 

In this paper we review the work of our groups on the 
application of the Vlasov kinetic equation to small- 
amplitude density vibrations in finite Fermi systems. 
The phenomenology that is described by this ap- 
proach includes giant resonances in nuclei, surface 
plasmons in atomic clusters and possible collective 
excitations of the electrons in heavy atoms (which 
however turned out not to be really collective). In 
the approach of Jeans , Vlasov |j and Landau 
iQ, ^ the theory has been developed for infinite homo- 
geneous systems, later Kirzhnitz et al. B extended 



it to non-uniform systems. 

The quantity that is studied in this approach is 
the phase-space density (or single-particle distribu- 
tion) /(r, p, i). This quantity, when multiplied by 
the phase space volume dvdp gives the number of 
particles contained in that volume and has several 
useful properties since its knowledge allows us to ob- 
tain information about macroscopic properties of the 
system like, for example, the ordinary density 



Q{r,t) = j dpf{r,p,t) 



or the pressure tensor 



Iljk{r,t) ^ J dppjVkf{r,p,t) 



(1) 



(2) 



We limit our interest to small oscillations of 
/(r,p,t) about its equilibrium value /o(r,p), which 
is supposed to describe the stationary state of the 
many-body system. If the system is subject to a weak 
external driving potential Vext{j^,t) — P{t)Q{r), then 
/o is changed by the small amount Sf{r,p,t), con- 
sequently the local equilibrium density also changes 

by 

6g{r,t)^ J dpSfir,p,t), (3) 
and the equilibriimr mean field Vo(r) by 



6V{r,t)= j dv'v{v,v')5g{Y',t). 



(4) 



The quantity -y(r, r') is the effective interaction be- 
tween two constituents of the many-body system. 
We shall consider only quantities at first order in 
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5f, moreover the equilibrium distribution will be as- 
sumed to depend only on the energy of the particles: 
/o(r,p) = F{ho{r,p)), where 



/io(r,p) = ^ + Vo{r) 
Zm 



(5) 



is the equilibrium hamiltonian. 

In a fully self-consistent approach the equilibrium 
mean field should be given by 

Voir) = J dr'v{r,r') J dp/o(r',p), (6) 

however we shall use phenomenological approxima- 
tions to it. There are two kinds of self-consistency 
requirements: the static self-consistency expressed 
by Eq. and the analogous dynamic self- 

consistency condition (^. We shall treat the static 
self-consistency condition more loosely, but will re- 
spect the dynamic condition because it is essential 
in the description of collective modes. 



2 HISTORICAL 
REMARKS 

The following differential equation for /(r,p,t) was 
considered by A. A. Vlasov in a study of plasma os- 
cillations 1^, ^: 



a/ 
dt 



m dr 



,dV. df 



dr dp 

together with the Poisson equation^ 



0. 



(7) 



(8) 



Because of Eq.(|l|), equations (|^) and (||) form a set of 
coupled equations that can be solved self-consistently. 
Thus the kinetic equation (|7|), supplemented by (^), 
can be considered as a dynamical generalization of 
the static Thomas-Fermi method. 



^Vlasov actually considered coupling to the full elecromag- 
netic field by using the Maxwell equations, here we simplify his 
argument and assume that the time dependence is sufficiently 
slow so that the laws of electrostatics can be applied at any 
instant t 



Althogh Eq.(0) is a simplified version of an equa- 
tion already derived by Boltzmann and the idea of 
self-consistency in connection with it had already 
been used by Jeans (see for a discussion of hys- 
torical priorities), we follow the use common both in 
plasma and in nuclear physics, and refer to Eq.(^ as 
"the Vlasov equation". 

The work of Vlasov was criticized by Landau |4j 
who pointed out some mathematical inconsistencies 
in Vlasov's solution and derived a rigorous solution 
of Eqs. for plasma oscillations. 

In his later work on Fermi liquids 1^ Landau ob- 
tained a kinetic equation similar to Eq.(|^) (see also 
Ref. 1^, p. 18) and used it to study the propagation 
of zero sound in liquid helium. 

Vlasov and Landau were interested in macroscopic 
systems, so they did not worry about surface effects; 
in their approach the system is supposed to be ho- 
mogeneous and infinite, so that translation invari- 
ance considerably simplifies calculations. Extension 
of their approach to microscopic systems like atoms 
or nuclei j^, apart from giving up translation in- 
variance, has to face the non-trivial problem of decid- 
ing which boundary conditions are most appropriate 
for the fluctuations of the single-particle distribution 
<5/(r,p,t). 

Kirzhnitz and collaborators |^ based their search 
for possible collective excitations of the electron cloud 
in heavy atoms on the linearized Vlasov equation. 
The main results of their work can be summarized 
in the two following equations for the polarization 
propagator: 

n(r,r',w) = n°(r,r',t^) (9) 
+ y dx y" dyn"(r,x,c^)z;(x,y)n(y,r',w). 



and 



n°(r,r',.) = -I^(.(r-r') (10) 



+iu; / dte-'"^' < 6{r{t) - r') > ) • 

The propagator n(r,r',a;) relates the (Fourier 
transformed in time) external disturbance at point 
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r' to the density response at point r: 



Sg{r,uj) 



(11) 



while n*' is the same propagator evaluated in the 
static mean field, that is, by neglecting the mean- 
field fluctuation induced by the external force. The 
most interesting part in the expression is the 
time integral. The quantity r{t) in the integrand is 
the coordinate of a particle that, at time t = is 
at point r and that moves in the static mean field 
according to the laws of classical physics. The S- 
function in the integrand contributes to the integral 
whenever r{t) = r', that is whenever a particle in r at 
t = has been through r' at some earlier time. The 
symbol < > in Eq. ( p^ ) means averaging over the 
directions of the particle velocity at t — 0, which re- 
quires that all the possible classical trajectories going 
from r' to r are taken into account when evaluating 
n". In Fig. 1 we give a schematic picture of a few 
classical orbits contributing to the time integral in 
Eq. ®. 

Bertsch |^ used the Vlasov equation as a start- 
ing point for a theory of giant resonances in nuclei. 
He pointed out that, since the Vlasov dynamics pre- 
serves the density of points in phase space, it does 
not violate the Pauli principle, even if it is a com- 
pletely classical theory. Thus, Eq.(0), after having 
been used by Jeans for stellar systems, by Vlasov 
and Landau for plasmas, by Landau for liquid he- 
lium and by Kirzhnitz et al. for atoms, found its way 
into nuclear physics as well. 

3 FIXED BOUNDARY 

3.1 Integrable systems 

The expression ( [To| ) for lP{r,r', uj) is quite general 
and in principle it is valid also for systems in which 
the motion of particles in the equilibrium mean field 
is chaotic, however it is not really suitable for prac- 
tical calculations since evaluating the time integral 
requires following the particle trajectories over an in- 
finite time interval. Kirzhnitz et al. avoided this dif- 
ficulty by restricting themselves to systems in which 






Figure 1: Schematic representation of a few classical or- 
bits contributing to the polarization propagator ( p^ ) for a 
finite system. Figure (a) shows the direct trajectory from 
r' to r, this is the only trajectory contributing to the 
propagator in a coUisionless homogeneous system. Fig- 
ures (b) and (c) show two typical trajectories with one 
and two reflections at the boundary. 
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the particle motion is periodic, but this is a very se- 
vere limitation. In the more recent work of Brink et 
al. it has been shown that there is an important 
class of finite systems for which the solution of the lin- 
earized Vlasov equation acquires a relatively simple 
form: they are all the systems for which the motion 
of particles in the equilibrium mean field can be de- 
scribed by an integrable Hamiltonian, this includes 
for example all spherically symmetric systems, but 
also some deformed systems. In these systems the 
particle motion is multiply periodic and, as a conse- 
quence, it is sufficient to follow the particle motion 
only over a finite time interval. The crucial step in 
applying the approach of to integrable systems, 
is an appropriate change of variables. The (r,p) 
variables are convenient for translation-invariant sys- 
tems, since in that case dVo/dr = and the lin- 
earized version of Eq. (|^) becomes simpler. A simi- 
lar simplification can be obtained for integrable sys- 
tems if action-angle variables (I, #) are used instead 
of (r, p). For these systems the action variables la 
are constants of the motion, while the conjugate an- 
gle variables are linear functions of time (see for 
example Ref. p. 457). An important property 

of these variables is that the motion is periodic in 
the angle variables with period 2tt. Consequently the 
field felt by a particle that is moving along such a 
trajectory can be Fourier expanded as 



(12) 



where n is a three-dimensional vector with integer 
components. Similar expansions hold also for the 
fluctuations of the mean field 

JT/(r,c^)=^5V„(I,c^)e^"-* (13) 

n 

and of the distribution function 

<5/(r,p,c.)=^<5/„(I,c.)e-*. (14) 

n 

When Eq.(^ is Fourier transformed in time and 
linearized by neglecting terms of order (Sf)"^ and 
higher, it reads 

-tL,6f{Lu) + {Sf,ho} (15) 
+F'{ho){ho, [V^'^'ioj] + SV{lo)]} = . 



Here the braces denote Poisson brackets, that can 
be evaluated according to any set of canonically con- 
jugate variables, using r and p gives the linearized 
version of Eq.(0), while using I and $ gives 



—iuj5f{uj) 
F'{ho)cj-( 



-u ■ 



Qyext Qgy 



(16) 



9* ) 



The three components of the vector u) determine the 
time-dependence of the angle variables: 



while, of course 



dia 



9feo(I.^) 



(17) 



(18) 



and hf) = /io(I). 

By using the expansions (12 1^), Eq.(|l^) gives 



5Ul,u) = F'{h„)\f3{Lo)Q^{l) (19) 
n • It) 



n • cj — (cj + ze) 



We have added a vanishingly small positive imaginary 
part ie in order to specify the behaviour of Sfn{to) at 
the pole f\ 

Equation (^9|) gives an explicit solution of the lin- 
earized Vlasov equation only if the fluctuation of the 
mean field SV can be neglected, otherwise it gives 
only an implicit solution, since 6V does depend on 
Sf; however its form immediately suggests an iter- 
ative procedure for obtaining the full solution: first 
evaluate Sf° by neglecting SV in (|l^, then evaluate 
SV by using Sf'^ as input and repeat until conver- 
gence is obtained. 

The quantity Sf is not very convenient for dis- 
cussing the solution of the linearized Vlasov equation, 
since it depends on the external driving force, so it 



^^In this subsection we have recalled the content of Sent. 3 n f 

fllCl]. Similar results had been obtained previously in Ref. |l3], 
A~D. apologises to those authors for not having been aware of 
their work before 
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is appropriate to introduce the polarization propaga- 
tors (^) and (|l^) , whose properties depend only on 
the system. Clearly, once we know Sf, by using Eqs. 
(|l|) and (pi]), we can obtain a corresponding expres- 
sion for the propagators II" and 11. We prefer to use 
their momentum representations, obtained from 
and ( p^ ) by taking Fourier transforms with respect 
to the coordinates: 



n(q,q',cj)= / dr / dr'e-*i '"n(r, r', w)e*i 



The expression of the 11° propagator is ul 



n°(q',q,c^ 
(2^) 



(20) 



(21) 



dIF'iho{l)) 
n-u){I) 



n • a;(I) - (lu + is) 
with the Fourier coefficients 



Q*M',l) Q„(q,I), 



Qn(q,I)=(^/d$e- 



(22) 



If the interaction i;(x,y) does depend only on the 
distance between particles, the integral equation (^ 
becomes Q 

n(q',q,c^) =n°(q',q,c^) (23) 
+ ^ / ^^knO(q',k,c^)i;(fc)n(k,q,c^). 

3.2 Spherical systems 

Spherical systems are an important class of integrable 
systems, hence they deserve a more detailed dis- 
cussion. Actually these systems are over-integrable, 
since in this case there are four constants of motion 
(the particle enegy e and the three components of its 
angular momentum A = r x p), so one of the an- 
gle variables must also be constant. The action- angle 
variables (I, $) allow for a very compact solution of 
the linearized Vlasov equation ( p^ ) , but they are not 
the only set of variables that can simplify the prob- 
lem. The important point is to choose as variables 
the largest possible number of constants of motion. 




Figure 2: Angular elements of the orbit. For a parti- 
cle moving in a central force field, the angles a and 
P are constant. 



so that the number of partial derivatives in Eq.(15) 
can be reduced. The variables used in instead 
of (r, p) are (e. A, r, a, /3, 7) of which four (e. A, and 
the two angles a and /3 shown in Fig. 2) are constants 
of motion. So we are left only with the two partial 
derivatives with respect to r and 7 in Eq.(p^. The 
7-derivative can be eliminated by means of the ex- 
pansion 



(24) 



LM 



and by using the well-known transformation property 
of the spherical harmonics Ylm{0, ^p) under the rota- 
tion specified by the Euler angles (a,/?, 7) (see e. g. 
Ill, p. 28): 



LM 



{9, v)^Yl (^MJv(«, 13, 7)) *Yln{0', v') . (25) 

N 



If the new reference frame rotates with the par- 
ticle, the angles 9' and i^' are constant and the 
only time-dependent angle is 7. The quantities 
2?£j^jY(a, /?, 7) are the rotation matrices and their ex- 
plicit 7-dependence, of the kind e~*^^, can be ex- 
ploited to eliminate the 7-derivative. In this way the 
initial three-dimensional problem is reduced to a one- 
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dimensional problem for particles moving in the ef- 
fective potential Vefj = Vo{r) + 

Then, as shown in the expansion ( p^ can be 
written as ^ 

Sf{r,p,u;)= (26) 



(l?^^(a,/3,7))*lL^(|,|), 



LMN 
X 



The functions Sf^^ refer to particles having positive 
('^/mat) '^^ negative (Sf^.j^) components of the radial 
velocity, of magnitude Vr{£, A, r) = ^/^(e — ^e//)- 

The linearized Vlasov equation ( [l5|) implies the fol- 
lowing system of coupled first-order differential equa- 
tions in the variable r for the functions Sf^j^ 



d_ 

dr 



B 



L± 
MN ■ 



with 
and 



iN 



X 



Wr(e, A, r) Vr{£, X,r) mr'^ 



-SMAr(e,A,r,w) 
d 



(27) 



(28) 



(29) 



iN 



A 



\or Vr(e, A, r) mr"^ 
l3{uj)QLM{r)+6VLM{r,oj) 



The functions QLRd{r) and dVLM{r,ui) are the coef- 
ficients of a multipole expansion similar to (|4|) for 
the external driving field and for the induced mean- 
field fluctuations, respectively. Because of Eq.M, the 
term 5VLA/(r, w) couples the two equations ( p7| ) for 
and Sf^^f^: 



SVLM{r,LLi) 



(30) 



8^2 



2L 



N=~L 



2' 2- 



■^Note that, unlike 



ike in |lo). 



here we do not include the factor 



de / dAA 



dr' 



TS^Lir, r') 



Vr{r') 

[SftMe, A, r', c.) + ^^/^-.^^(e, A, r', c^)] . 

Before solving the system of equations (p7|), we 
must specify the boundary conditions satisfied by 
the phase-space density fluctuations. The conditions 
used in [flO| are 



SfMNin) 



'^/mw('^2), 



(31) 
(32) 



where ^1(2) are the classical turning points, for which 
Wr(e, A, r) = 0. 

With these boundary conditions the solution of 
the linearized Vlasov equation for spherical systems 
agrees with that given by the method of action-angle 
variables [|lOj. 

There are some some further simplifications for the 
polarization propagators of spherical systems, com- 
pared to those given in the previous subsection. The 
three-dimensional integral equation (|2^) reduces to 
a set of one-dimensional integral equations for each 
multipolarity L 



IlL{q',q,u;) = nl{q',q,Lu) + 



1 



(27r)3 

dkk^ nl{q\ k, uj) v{k) WL{k, q, u) , 



(33) 



and the 11° propagator also becomes somewhat sim- 
pler than (|l]): 



nl{q',q,Lu)^ 

+ 00 



(34) 



2L + I E E r^^(2 2- 

n— — 00 N— — L 

dcF'ie) fdXXT- ^^o + Nu;, 
J ? 



nujQ + Nuo^p — (ix) + ie) 



Q[^;>\q\e,X)Q'^j>{q,e,X). 

Here LOQ[t,X) is the frequency of radial motion of 
a particle with enegy e and magnitude of angular 
momentum A in the effective potential Ve//, while 
a;^(e. A) is the precession frequency of the periapsis 
in the plane of the orbit (jl^,p.509). These two fre- 
quencies determine the shape of the orbit in a central 
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potential. The period of radial motion T is T(e, A) = 
2tt /ijjQ. The Fom'ier coefRcients Ql^j^{q, e, A) are given 

by 

Qi^^(9,e,A)= (35) 

2 dr 

/ s cos[s„Af(r)]Q^(gr), 

with 

gL(gr) =4^z^jL(qr). (36) 
The phases SnN{f) are 

s„Ar(r) = nuiQTir) + N{uj^T{r) - 7(r)) , (37) 

where 

, . r dr' 

is the time taken by a particle to move from the inner 
turning point ri to point r, and 

7(^)=r^^A (39) 
Jj._^ Vr[e, A, r'j mr'^ 

is the angle spanned by the radius vector during the 
time interval T(r). 

In Eqs. (H^) and ( p4[ ) it has been asumed that 
the equilibrium distribution of particles depends only 
on the particle energy. Usually one takes F{e) oc 
0(ei?—e), where ep is the Fermi energy, so that F'{e) oc 
5{ep — e) and the expression ( |3^ ) becomes simpler. 

At first sight it might seem that the infinite sum 
over n in ( ^4| ) might create difficulties, but in practice 
it is enough to include in the sum the very first few 
terms around rt = to get a sufficient approximation. 

3.3 Atomic plasmons? 

The early calculations of Kirzhnitz and collabora- 
tors suggested the existence of two collective high- 
energy excitations in the spectrum of heavy atoms 
1^ . These excitations had been interpreted as collec- 
tive, plasmon-type oscillations of the atomic electron 
cloud. In order to apply their formalism, Kirzhnitz 
et al. had to approximate the atomic self-consistent 
mean field with a potential in which electrons hav- 
ing zero total energy move along closed trajectories. 



Our closely related approach allows also for orbits 
that are not closed, but, of course, it can be applied 
also to closed trajectories. This has been done in 
p^ , where the calculations of Kirzhnitz et al. have 
been repeated by using the formalism discussed in 
this Section. In that study only one "collective" solu- 
tion of the integral equation ( |3^ ) was found, moreover, 
thanks to the analytical insight given by the semiclas- 
sical method, this solution was shown to be qualita- 
tively different from the plasma oscillations of a uni- 
form electron gas. While plasma oscillations belong 
to a strong-coupling regime (see e.g. @, p. 185), 
the "collective" solution found in pertains to a 
regime of weak coupling. Thus it was concluded that 
there is no evidence for plasmon-like excitations in 
this semiclassical theory of atomic spectra. This con- 
clusion is in agreement with the result of analogous 
quantum calculations. 

3.4 Integrable spheroidal systems 

The equilibrium self-consistent mean field in heavy 
closed-shell nuclei is well approximated by a spher- 
ical Saxon- Woods potential, with parameters that 
can be determined phenomenologically. A somewhat 
similar mean field describes also metal clusters |l8| . 
Thus, for open-shell systems it is reasonable to try a 
deformed Saxon- Woods like potential. However the 
results of the first subsection cannot be applied to 
this kind of potentials, since accurate calculations of 
the classical phase space for axially deformed Saxon- 
Woods potentials have shown that it contains regions 
of chaoticity ||l^ . But, if the surface diffuseness is ne- 
glected and a spheroidal cavity is assumed to approx- 
imate the mean field of large deformed nuclei or 
clusters, then that formalism can be applied directly. 
This has been done in Ref. to study the surface 
plasmons in deformed atomic clusters. The main re- 
sult of that paper is that, under certain conditions, 
oblate deformed clusters may display a more complex 
shape of the surface plasmon peaks than predicted by 
the classical Mie theory (three peaks instead of two) . 
This prediction has not yet been checked experimen- 
tally 

Fig. 3 shows essentially the single-particle (shaded 
black) and collective (shaded grey) dipole strength 
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Prolate Cavity 




Figure 3: Photoabsorption cross section per valence elec- 
tron (shaded grey) for prolate and oblate sodium clus- 
ters. The black-shaded peaks show results obtained in 
the single-particle approximation in which the mean-field 
fluctuation ^ is neglected. The arrows indicate the po- 
sition of the plasmon peak in the spherical case. 



functions for prolate and oblate spheroidal clusters. 
The figure is taken from Ref. The deforma- 

tion parameter r/ is the ratio between the larger and 
smaller radius of the spheroid. The frequency is ex- 
pressed in units of the Mie frequency lomib- In the 
lower panel the fragmentation of the high-frequency 
peak is clearly visible for deformations rj =^ 1.25 and 
T] = 1.5. 

4 MOVING 
BOUNDARY 

4.1 Connection with liquid drop 
model 

In heavy closed-shell nuclei the equilibrium mean 
field can be further approximated by a spherical 
square-well potential of radius R « 1.2A3fm. In this 
simplified picture several of the quantities needed to 
evaluate the linear response of the system can be ex- 
pressed analytically, hence it gives a useful reference 
model. 

When the surface diffusion is neglected, the bound- 
ary condition (^) corresponds to requiring that the 
particles collide elastically with the static nuclear sur- 
face r2 — R (mirror-reflection condition). 

In Ref. it has been shown that the results of 
pod can be generalized to a model in which the sharp 
surface is allowed to vibrate, simply by changing the 
boundary condition (^). 

If the radius of a spherical nucleus is allowed to 
vibrate according to the usual liquid-drop model ex- 
pression Q 

R{9,ip,t) = R + Y, ^Rlm {t)YLM iO,ip), (40) 

LM 

then the mirror-reflection condition (|3^) is changed 
to 

SftmiR) - SfNmiR) - 2F'{e)iupr5RLM[u;) , (41) 

where Pr is the magnitude of the radial component 
of the particle momentum. We put a tilde over the 
moving-surface solutions to distinguish them from 
the untilded fixed-surface solutions. The boundary 
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condition ( [4lD still corresponds to a mirror reflection, 
but in the reference frame of the moving surface. 

It can be easily checked by direct substitution that 
the following integral equation 



with 



andQ 



(42) 



F'(e 



,±i[u>T{r)-N'l{r)] 



-UjprSRLM{uj) 



sm[ujT{R) - N-f{R)] 

Jri 



A 



-SA/Ar(e,A,r,w) = 

\or Vr(e, A, r) mr"^ 
l3{uj)QLM{r)+6VLM{r,uj 



6VLM{r,uj) ■ 

.2 L 



2L 



f f dr' 
de / dAA / — rT:'"L{r,r') 

(e, A, r',iu) + 6f^fj,{e, A, r', cu)] 



|gi2[c^r(i?,)-Af7(fl)] 

is equivalent to the system 



L± 

MN 



(43) 



(44) 



(45) 



(46) 



^There is a misprint in Eq. (3.16) of Rcf . [psj , the present 
expression is the correct one. 



with the boundary conditions ( |3l| ) and (^), and that 
it respects the self-consistency condition (|). 

Equations ( ^^ - |4^ ) give an implicit solution of 
the linearized Vlasov equation with moving-surface 
boundary conditions. They are more complicated 
than the corresponding fixed-surface equations be- 
cause they contain also the unknown collective co- 
ordinates SRlm{^)- In order to proceed further, we 
need an additional relation that specifies these quan- 
tities. As shown in this relation can be found 
within the liquid-drop model if one recalls that, in 
that model, a change in the curvature radius of the 
surface results in a change of the pressure at the sur- 
face given by (cf. Eq. (6A-57) of H) 



6V{RA^,^)-T.cJ-^^Y,m( 

LM 



,ip). (47) 



The restoring force parameter Cl can be easily 
evaluated: if the Coulomb repulsion between protons 
is neglected. 



Cl = cri?2(L- l)(L-|-2), 



(48) 



where a « IMeVfm"^ is the phenomenological 
surface-tension parameter, while taking into account 
also the Coulomb interaction changes somewhat this 
expression Q. 

The connection with kinetic theory is then made 
simply by noticing that the surface pressure (B^) is an 
appropriate diagonal element of the tensor (0) (gen- 
eralized to Fermi liquids, see e. g. p4]): 



SV{R,0,(p,Lu) = lim 

r—^R 



(49) 



dpprVr{Sfir,p,u;) - F'{e)SV{r,u;)) 



This integral must be evaluated for r < R and then 
we must let r — > i? from inside. 

Multiplying both (|^ and (|^) by y£^,f(6',(^) and 
integrating over the directions of r, gives 



SRlmIi^) 
87r2 i?2 



(50) 



StT R sr-^ I '''x 
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de / dXXpr 



with 



~2F'{e)SVLM{R,uj) 
with the definition 



SftiN^Sf^+, + Sf^f^. (51) 

Equation ( |50| ) relates the fluctuations of surface and 
distribution- function in the approach of ; it is an 
imphcit relation for SR, since both Sf and SV do 
depend on SR. 

4.2 Approximate solutions for surface 
vibrations 

In order to obtain an explicit expression for the col- 
lective coordinates 5Rlm{'^) we make the following 
approximation: 



(52) 



-2F'{e) 



cos[aJT(r) — N^{r)] 
(R) - Nj{R)] 



sm CiJT 



UjprSRLM{i^) , 



with 6f^jj^{e, X,r,uj) the fixed-surface solution. This 
approximation is based on the assumption that, apart 
from the term proportional to SR^m, which alone 
fulfills the boundary condition (|4l|), in the bulk of 
the system the moving-surface solution does not differ 
too much from the fixed-surface one. 
Then, Eq. (Q) gives: 



5Rlm{^) 
8^2 



(53) 



{ 



L 



2L 
de 



1 ^ |^^^(2'2^ 



dXXpr 

-2F\e)SVLM{R,Lo)\}[^ 

t |n»(f,|)i7*n.) 



N=-L 



dXXpr ^ujpr cot[a;r(i?) — N^{R)] 

iR,^)]Y\ 



(4^)^ 

n, 

2L + 1 



(54) 



N=-L 



2' 2- 

deF'{e) J dXX J dr'vL{r,r') 

cos[a;r(r') — N'y{r')] 
sin[wT(i?) - Nj{R)] ^ ■ 

Equation ( ^3| ) is an interesting, although approxi- 
mate, expression for the collective coordinates given 
by the approach of Ref. j2^ . The vanishing of the 
denominator in this expression determines the poles 
of (5i?LM(w), hence the eigenfrequencies of the sur- 
face vibrations 6RLM{t)- Note that, according to Eq. 
([5^), for L 7^ 1 these eigenfrequencies depend on the 
surface tension, which is physically sound. 

The eigenfrequencies determined by the vanishing 
of the denominator in Eq. ( p3| ) can be directly com- 
pared to the experimental excitation energies of the 
corresponding modes. Up to now this has been done 
only after making an additional approximation: if 
the fixed-surface solution is approximated by its zero- 
order expression Sf1/"j^ , appropriate for a gas of non- 
interacting particles, a corresponding zero-order ap- 
proximation is obtained for SR^Mi^) 



srIaA^) = 

N=-L 

de I dXXpr[df'i/^j^ie,X,R,iu) 



(55) 



2L + 



7\7 T ^ 



N=-L 



dXXpl [ coi[uJT{R) - iV7(i?)] | 



Note that this approximation contains more dynam- 
ics than the corresponding fixed-surface zero-order 
approximation, since the interaction between parti- 
cles is taken into account to some extent through the 
phenomenological surface-tension parameter a. 
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4.3 Surface and compression modes 

It is interesting to compare the dynamics described 
by Eq. ( |55| ) with that of the liquid-drop model. 

The vanishing of the denominator in Eq. (55) can 
be written as 



-0, 



(56) 



i?2 i?2 

with the function Xl{^) defined as in Ref. [ p5[ : 

XlH = (57) 

,2(4vr)^ 
1 

N=-L 



ujR 



2L 

N=-L 

deF'ie) [ dWpl cot[u;T{R) ~ N-f{R)] . 



The expression analogous to (pq ) in the liquid-drop 
model is 

Cl - lo^Dl = , (58) 

or, allowing also for compression modes (cf. Eq. (6A- 
58) of H), 



0. 



(59) 



\r=R 



with Uc the velocity of sound. This condition deter- 
mines the eigenfrequencies of the various modes in 
the liquid-drop model |Q. These eigenfrequencies 
depend also on the mass parameters Dl that can be 
evaluated within that model, but one has to make 
some assumption about the kind of flow occurring in 
the system p3{ . 

In order to make a comparison between the solu- 
tions of (^) and ( |59| ) it would be desirable to have a 
simple analytical expression for the function Xi(w). 
The pole expansion of the cotangent 



cot z 



1 



71— — OO 



(60) 



allows us to express Xl{^) in £^ form that is more 
similar to the 11^ propagator (^3): 



XL(t^) = ujR 



2 (iH 

2L + 1 



(61) 



OO L or 

n=-ooN=-L 

dXXpl - 



' {uo + ie) - {nujo + Nui^) ' 



but this only shows that XL(t^) is a rather involved 
complex function. The very fact that xl is complex 
is already an interesting result since we can expect 
complex eigenfrequencies as solutions of Eq. ( |56| ) , in 
close analogy with the Landau damping phenomenon 
in homogeneous system and in contrast to the liquid- 
drop model condition ( |59| ) that involves only real 
quantities. 

In the case L ~ the function Xi(^) can be evalu- 
ated analytically while for L = 1,2,3 relatively 
simple expressions have been obtained for the coeffi- 
cients of the low-frequency expansion 



(62) 



For isoscalar monopole vibrations, the solution of 
Eq.(|5^) compares reasonably well with the experi- 
mental value of the monopole giant resonance energy 
in heavy nuclei. 

For isoscalar quadrupole oscillations, in it was 
concluded that the approximation ( p5| ) is not ade- 
quate to describe the low-energy isoscalar quadrupole 
excitations. It is not clear if this failure is due to the 
neglect of the bulk mean-field fluctuations in (^) or 
to some other reason. Further work is needed on this 
problem. 

About the low-energy isoscalar octupole oscilla- 
tions, in 1 25 it was concluded that they are over- 
damped in the approximation given by Eq. (^5|), 
and this might offer a qualitative explanation for the 
background observed in inelastic proton scattering 

il- 

Another application of Eq. (55) has been made 
in Ref. |2^. There this equation has been used to 
study the isoscalar dipole compression mode. In or- 
der to get reasonable results in this case it is vital to 
treat the centre of mass motion correctly, and this can 
be done in a fully self-consistent approach. Because 
Eq.(|5^) is not really self-consistent, in |^ it was 
found necessary to add an extra term to the dipole 
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10 20 30 

E (MeV) 



Figure 4: Intrinsic dipole energy-weighted strength 
function compared to experimental data for ^°^Pb 
(a), "^Sn (b) and 9°Zr (c). The soHd curve is for a 
confined gas of A non-interacting nucleons with in- 
compressibility K « 200MeV, the dashed curve is for 
interacting nucleons with K = 160MeV. The data 



are from Ref.|28 



response function based on Eq. (pq). For other mul- 
tipolarities this extra term can also be interpreted as 
the contribution to the response due to the change of 
shape of the system in the moving-surface approach, 
for L = 1 the system does not change its shape, but 
can translate as a whole. Since there is no force op- 
posing this translation, this centre of mass motion 
results in a pole at uj — in the dipole response 
function. The interesting physics is contained in the 
intrinsic response, that is shown in Fig. 4. 

The results of the moving-surface theory are in 
qualitative agreement with experimental data on the 



isoscalar dipole mode, as can be seen from Fig. 4. 
Moreover, the insight given by this semiclassical ap- 
proach gives us the possibility to establish a direct 
link between the peak profile of the isoscalar dipole 
resonance and the incommpressibility of nuclear mat- 
ter ||2^. This link is evidenced by the dashed curve 
in Fig. 4. 

The approach described here can be easily general- 
ized to a two-component fluid. In p9[ | such a gener- 
alization has been used to study the effect of neutron 
excess on the isoscalar and isovector monopole re- 
sponse functions. For such a two-component system, 
Eq.(^) is replaced by the two following equations 



and 



n;v(i?,f?,^,L^)) 



(63) 



(64) 



Here Ilrr^'' (i?, 9, if, u) is the normal component of the 
momentum-flux tensor [r.h.s. of Eq. (|49|)] associated 
with neutrons (protons) only. Apart from the usual 
surface pressure 6V that is present when the neutron 
and proton surfaces move in phase, in this case there 
is also an additional pressure SVq that is caused by 
the forces opposing the pulling apart of the neutron 
and proton surfaces when they move out of phase. 
This extra pressure can also be related to an appro- 
priate phenomenological parameter. With the sur- 
face energy of [|o|, the pressure 5Pq can be written 
as 



SrQ{R,e,^,u;) = 

-^[SRni9,ip,iv) -6Rpi9,ip,iv)] , 
ztttq 



(65) 



where Q is the neutron skin stiffness coefficient of 
that is analogous to the surface tension parameter a 
of Eq. (g|) and tq = 1.14fm. 

As shown in j29j, the neutron excess does not af- 
fect much the isoscalar giant monopole resonance, 
while its effect on the isovector resonance is more 
pronounced. 
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Figure 5: Isoscalar dipole strength functions eval- 
uated in zero-order approximation with moving- 
surface (fuU curve) and fixed-surface (dashed curve) 
boundary conditions. The fuU curve corresponds to 
the full curve in Fig. 4 

4.4 Fixed vs. moving boundary 

The possibility of solving the Vlasov equation with 
different boundary conditions allows us to decide 
which condition is more appropriate by comparison 
with experiment. In Fig. 5 we give an example of a 
nuclear response function evaluated according to the 
fixed- and moving-surface boundary conditions. 

In this case there is no doubt that comparison with 
data (shown in Fig. 4) favours the moving-boundary 
solution, at least at the level of zero-order approx- 
imations. The present example is a rather extreme 
case, for other response functions the difference be- 
tween the two solutions is not expected to be so large. 
Work on detailed comparison between the two solu- 
tions for different multipolarities is still in progress. 



5 QUANTUM EFFECTS 

The Vlasov theory is completely classical, however it 
is easy to establish a connection with quantum me- 
chanics. It is well known that the time-dependent 
Hartree-Fock equation for the Wigner transform of 
the quantal one-body density matrix reads (see e.g. 
Ref. m, p. 553) 



1 

X sm[-h{V -^p-^p- ^)] 
+F(r,p,t)] =0 



(66) 



and reduces to 



dt ni dr 

,dV, df ,dV, df „ 
-( — ) • — + ( — ) • — = 
dr dp dp dr 



(67) 



if ?i ^ 0. The Hartree-Fock potential V{r, p, t) 
is momentum-dependent because the Fock term is 
non-local in coordinate space. If the momentum- 
dependence of V is neglected (Hartree approxima- 
tion), then Eq. (^ coincides with Eq. (^. 

There are several possibilities for including quan- 
tum effects into the Vlasov theory. The most trivial 
one is to choose an appropriate equilibrium distri- 
bution F{ho{r,p)). The standard choice for zero- 
temperature fermions is 



Fiho{r,p)) 
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{2'Kny 



:'d{eF - ^o(r,p)) : 



(68) 



where ei? is the Fermi enegy and g is the degeneracy 
factor (g = 2 (4) for electrons (nucleons)). 

The choice (^8|) respects the Pauli principle in the 
sense that there cannot be more than one fermion 
of given spin (and isospin) in a phase-space cell of 
volume . As stressed by Bertsch Q, the Liouville 
theorem then ensures that the Pauli principle will be 
respected during the time evolution of the system de- 
scribed by the Vlasov equation. Thus the Liouville 
theorem, stating that the points representing the sys- 
tem in the 6-dimensional phase space evolve as an 
incompressible fluid, justifies the use of the Vlasov 
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equation for quantum systems. Using the equilib- 
rium distribution ( p8| ) makes the completey classical 
Vlasov approach a semiclassical one. Clearly in this 
way we do not expect to reproduce details, but only 
gross properties of quantum many-body systems. 

In principle one could introduce further quantum 
effects also by taking into account terms of order 
and higher that have been neglected in Eq. (|67|). The 
problem with such a direct approach is that the semi- 
classical calculations become quickly more compli- 
cated than the corresponding quantum calculations, 
so this method would not be very convenient. A more 
pragmatic attitude has been taken in Refs. and 
p3[ , where the semiclassical propagators given by the 
Vlasov theory have been compared directly with the 
analogous quantum propagators. 

In the relation between the Fourier coefRcients 
( p^ and the corresponding quantum matrix elements 
has been discussed. If single-particle matrix elements 
are evaluated in the WKB approximation, they tend 
to Fourier coefficients analogous to ( p^ ) in the limit 
of large quantum numbers, but this is not necessarily 
the case for the exact quantum matrix elements. 

In the closed-orbit theory, that had been used 
by various authors to evaluate quantum corrections 
to the Thomas-Fermi level density, has been applied 
to evalute quantum corrections to the Vlasov prop- 
agator (H^). The resulting strength function com- 
pares better with the analogous quantum strength 
function. A fascinating aspect of this study is the 
possibility of relating features of the excitation spec- 
trum of quantum systems to simple properties of the 
classical orbits. 

6 CONCLUSIONS 

In the kinetic-theory approach to linear response dis- 
cussed here the finite size of the system has been 
explicitly taken into account. Two different bound- 
ary conditions for the density oscillations have been 
considered: fixed and moving surface. 

The fixed-surface solution is closely related to the 
quantum random phase approximation (RPA) since 
the integral equation for the semiclassical collec- 
tive propagator is essentially the same as the equa- 



tion for the quantum RPA propagator when exchange 
terms are neglected (Hartree approximation). The 
main difference between the semiclassical and quan- 
tum approaches is in the "free" propagator 11°. While 
the quantum expression for this propagator involves 
wave functions and single-particle energies, the semi- 
classical propagator involves only properties of the 
classical orbits and is considerably simpler to evalu- 
ate. 

The moving-surface solution has been related to 
the liquid-drop model of nuclei. Compared to that 
model, the present approach is more microscopic and 
it allows also for the possibility of Landau damping 
of collective vibrations. 

When compared to experiment, in some cases the 
semiclassical solutions have the same problems of the 
RPA and of the liquid-drop model. For example, the 
fixed-surface approach is, like RPA, able to describe 
reasonably well the position of the surface plasmon 
in atomic clusters, but, again like RPA, fails to re- 
produce the observed width of this resonance. Possi- 
ble moving-surface contributions to this width have 
not yet been investigated in detail within the present 
approach. The moving-surface approach shares with 
the liquid-drop model the difficulties in reproducing 
the experimental value of low-energy quadrupole os- 
cillations in heavy nuclei. However the insight given 
by this approach could be useful for improving both 
the quantum RPA theory and the liquid-drop model. 
In conclusion, the semiclassical theory discussed here 
offers an interesting alternative to the quantum RPA 
approach because it is simpler to implement numer- 
ically, but it is still sufhciently sophisticated to take 
into account some realistic features (finite size, equi- 
librium shape) of the many-body system. 
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